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We present new Neutron Spin Echo (NSE) results and a revisited analysis of historical data 
on spin glasses, which reveal a pure power-law time decay of the spin autocorrelation function 
s(Q,t) — S(Q,t)/S(Q) at the glass temperature T 9 , each power law exponent being in excellent 
agreement with that calculated from dynamic and static critical exponents deduced from macro- 
scopic susceptibility measurements made on a quite different time scale. It is the first time that 
this scaling relation involving exponents of different physical quantities determined by completely 
independent experimental methods is stringently verified experimentally in a spin glass. As spin 
glasses are a subgroup of the vast family of glassy systems also comprising structural glasses and 
other non-crystalline systems the observed strict critical scaling behaviour is important. Above the 
phase transition the strikingly non-exponential relaxation, best fitted by the Ogielski (power-law 
times stretched exponential) function, appears as an intrinsic, homogeneous feature of spin glasses. 
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The glass transition, characterised by a dramatic slow- 
ing down of the dynamics without any noticeable change 
in the spatial order, is a generic phenomenon, seen in sys- 
tems as different as disordered magnets, polymers and 
biological substances. In spite of its universality and 
of intense experimental and theoretical efforts it is still 
controversial whether the glass transition in structural 
glasses is a gradual freezing or a phase transition. The 
difficulty to identify the nature of the glass transition 
is due to the absence of an observable order parame- 
ter analagous to magnetization in the low temperature 
phase, usually a key quantity in the study of phase tran- 
sitions. This is due to the absence of any static spa- 
cial fingerprint; instead, the order parameter appears in 
the dynamics^. In fact, the "snapshot" structure factor 
S(Q) = S(Q, t = 0), which reflects the short and medium 
range static correlations, shows no essential change when 
passing from the high temperature liquid (or paramag- 
net) to the low temperature frozen glass phase. In this 
situation the observation of dynamic scaling relations, 
which are the direct consequence of the homogeneity hy- 
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pothesis in the vicinity of a critical instability 3 can reveal 
the crucial signature of a true phase transition. 

In spin glasses, which are the simplest realisations of 
glassy systems from the experimental as well as from 
the theoretical point of view, a phase transition is well 
established^. Low frequency dynamic susceptibility as 
well as the non linear part of the static susceptibility fol- 
low the scaling relations and the analysis leads to accu- 
rate determination of the static critical exponents 7, /3, 8 
and of the dynamic exponent a&LSiS. The verification of 
scaling relations between certain exponents, determined 
by completely independent experimental methods, how- 
ever, was up to now impossible due to the absence of 
any obvious critical behaviour on other physical quan- 
tities, like the specific heat. Here we show that scal- 
ing relations can also be verified experimentally in spin 
glasses. New Neutron Spin Echo (NSE) results and also 
a revisited analysis of historical datai2*ii*i2i show a pure 
power-law time decay of the spin autocorrelation function 
s(Q,t) — S(Q,t)/S(Q,t = 0) at the glass temperature 
T g , with a power law decay exponent which is in excel- 
lent agreement with that calculated from the dynamic 
and static critical exponents deduced from zero and low 
frequency susceptibility measurements. The interplay 
between neutron scattering, macroscopic magnetic, hy- 
perfine field measurements, and simulations has always 



been decisive in understanding spin glasses and these re- 
sults constitute the strongest experimental evidence yet 
for a true phase transition with a non-conventional order 
parameter; they also imply that the prominently non- 
exponential relaxation is an intrinsic, homogeneous fea- 
ture of spin glasses. 
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FIG. 1: Temperature dependence of the normalised interme- 
diate scattering function s(Q,t) of Auo.86Feo.14. The spec- 
tra were collected at Q = 0.4 nm _1 with the Neutron Spin 
Echo spectrometer IN15 (ILL) for T=30.7 K (closed squares), 
40.6 K (~T g , open squares), 45.7 K (closed circles), 50.8 K 
(open circles) and 55.8 K (closed rhombs) respectively. The 
continuous lines are the best fits to the data of a simple power 
law decay below T 9 (~41 K) and of the Ogielski function 
above T s (see text) 

Most of the information about glass transitions comes 
from the high temperature unfrozen phase, where ther- 
modynamic equilibrium is easily reached without any 
long-time drifts and aging phenomena. In spin glasses, 
where the magnetization is always zero, the fundamen- 
tal parameter is the mean spin autocorrelation function 
q(t — t ) =< Sj(t).Sj(t ) > where Si is the spin at a site 
i and the average runs over all sites and configurations 
of the sample. Critical behaviour in the paramagnetic 
phase is seen in the non-linear susceptibility (or "spin 
glass susceptibility" ) Below T g the Edwards Ander- 

son order parameter q(t— s-oo) = <Sj(t = 0)Si(t— >oo)>, 
becomes non-zeroii^ Neutron Spin Echo spectroscopy 
measures the scattering function S(Q,t) and after nor- 
malization by S(Q,0) delivers a direct determination of 
the autocorrelation function q(t) = s(Q— >0,£). NSE cov- 
ers a time domain ranging from 10 -12 s to some 10~ 8 s, 
i.e. from characteristic microscopic times up to times, 
which already belong to the "long" time relaxation do- 
main. The first NSE experiment ever performed on a 
glassy system was made on the reference spin glass CuMn 
5% in 1979iS and the results strongly influenced subse- 
quent thinking on (spin)glass dynamics. For the first 
time is was shown that non-conventional dynamics is not 
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FIG. 2: NSE spectra of Au . 8 6Fe .i4 collected at Q=0.08A _1 
at the ILL spectrometer IN15 (full symbols) and at the 
BENSC spectrometer SPAN (open symbols) plotted in a lin- 
log scale. The circles are measured at 45.6 K and the triangles 
at 40.6 K respectively. For the sake of clarity the data sets 
corresponding to each of the temperatures have been shifted 
with respect to each other in the vertical scale. The data 
at 40.6 K are fitted to a simple power law (continuous line). 
The continuous line through the 45.6 K data represents the 
fit of an Ogielski function (see text) . The dashed and dotted 
curves correspond to a stretched exponential and a simple 
exponential decay respectively. 



limited to the spin glass phase but also extends into the 
paramagnetic phase well above T g . Non-exponential and 
Q-independent relaxation occurs in a large temperature 
range up to 2-3 T g , which can arguably be identified with 
the Griffiths phased. For about T>1.2 T g the relaxation 
can be described by a broad distribution of Arrhcnius 
activation energies. Closer to T g , however, a more dra- 
matic slowing down sets in, which can be interpreted as 
the footprint of a phase transition with a critical region of 
usual extent. Here we report on the first detailed analysis 
of s(Q, t) around T g in spin glasses, based on enhanced 
quality data obtained by using new generation NSE spec- 
trometers. 

For an accurate determination of the NSE spectra we 
chose Auo.86F e o.i4- AuFe is a classical metallic Heisen- 
berg spin glass with significant local anisotropyi^ and 
with strong ferromagnetic correlations which amplify 
the magnetic scattering in the forward direction so im- 
proving the ratio between the magnetic signal and all 
non-magnetic (structural) contributions, i.e. the signal 
to noise ratio. The sample was a polycrystalline disc 
0.5 mm thick with a diameter of 37 mm prepared by 
arc melting of the constituents. It was subsequently 
cold worked, homogenised at 900°C, annealed at 550°C 
and then quenched and kept in liquid Nitrogen?lS Given 
the vicinity to the percolation threshold x c (~15.5% at. 
Fe), above which ferromagnetism sets in^£ the anneal- 
ing and quenching procedure was repeated before every 
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FIG. 3: s(Q,t) of Auo.86Feo.i4 measured at 0.04 A -1 (close 
symbols) and 0.08 A" 1 (open symbols) for 45.6 and 40.6 K 
respectively. For the sake of clarity the data sets correspond- 
ing to each of the temperatures have been shifted with respect 
to each other in the vertical scale. 



series of measurements. The NSE data were collected at 
the high resolution spectrometer IN15 of ILLi& at an in- 
coming wavelength of 0.8 nm for Q = 0.4 and 0.8 nm -1 
respectively. These results were supplemented by mea- 
surements at the wide angle NSE spectrometer SPAN 
of BENSC— at an incoming wavelength of 0.45 nm for 
0.6 nm _1 <Q < 2.6 nm -1 . We used the paramagnetic 
NSE set-up, which directly delivers the magnetic part of 
the NSE signal and for this reason no background cor- 
rection was required. All NSE spectra were normalized 
against the resolution function of the spectrometers, de- 
termined with the sample well below T y , at 2 K, where 
the spin dynamics is completely frozen. A small part of 
the sample was taken out for dc susceptibility measure- 
ments with a commercial SQUID magnetometer at the 
HMI and for ac susceptibility measurements, which were 
made from 10 Hz up to 10 KHz with the MAGLAB setup 
at the Physics and Astronomy Department, University of 
Leeds. The spin glass temperature of T g = 41.0 ± 0.3 K 
was determined from the maximum of the static suscep- 
tibility. 

The normalised intermediate scattering function s(Q,t) 
of Auo.86Feo.14 at Q = 0.4 nm -1 is shown in Fig. 1 
plotted in a log-log scale. The spectra span a dynamic 
range of three orders of magnitude and by combining 
spectra collected at two wavelengths on INI 5 and SPAN 
the time domain of the observation is extended up to 
almost 4 decades (Fig. 2). The time dependence of the 
experimental s(Q,t) is impressively similar to that of the 
numerical q(t) found in large scale Ising spin glass simu- 
lations, which revealed the existence of a phase transition 
in three-dimensional Ising spin glasses. 20 From quite gen- 
eral scaling argumentsj2i at a continuous phase transition 
relaxation must be of the form t~ x f(t/r(T)) where r(T) 
diverges as (T — T y )~ z , and / is a non- universal func- 



tion to be determined for each system. For a spin glass 
the power law exponent x is related to the standard crit- 
ical exponents through x = (d — 2 + r;)/2a2&. Here r\ 
is the Fisher or "anomalous dimension" exponent and z 
is again the dynamical exponent. The Ising simulations 
showed that, as T g is approached from above, q(t) is 
strongly non-exponential. Ogielski chose to represent / 
by the stretched exponential or KWW function, familiar 
in fragile glass dynamics. Excellent fits were obtained 
with q(t) cx t~ x exp((—t /r(T)) /3 ) and T dependent r and 
(3. (The KWW (3 is not to be confused with the crit- 
ical exponent (3). He and others found a temperature 
dependent (3 tending to near 1/3 at T g and increasing 
with Tm^^^ The most important point for our data 
analysis is that precisely at T g dynamic scaling predicts 
a pure power law decay for the autocorrelation function : 
q(t) cx t~ x . This rule is quite general; its functional form 
does not depend on details such as the Ising character of 
the spins but the value of x depends on the exponents 
for the particular system under study. Scaling therefore 
gives a unique opportunity for describing the NSE spec- 
tra at T g in terms of exponents determined by completely 
independent dynamic and non-linear macroscopic suscep- 
tibility measurements. 
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FIG. 4: Combined NSE spectra with s(Q=0,t) values deduced 
from macroscopic ac susceptibility measurements below T 9 . 
Also in the close vicinity of T 9 , the pure power law decay of 
s(Q,t) holds over an impressively large dynamic range of more 
than nine orders of magnitude in time, from the microscopic 
to the macroscopic times. 

As seen on Fig. 1, for T<T S we found s(Q,t) cx t~ x 
and the lines represent the best fits to the data with 
x=0.116±0.007 and 0.025±0.005 at T=40.6(~T 5 ) and 
30.7 K respectively. Above T g the relaxation is also 
strongly non-exponential. As shown in Figure 2, pure ex- 
ponential as well as stretched exponential decay (without 
the power-law prefactor) can be definitively ruled out at 
all T in the range studied. In fact at 45.6 K, the fit with 
the Ogielski function leads to \ 2 — 0.28, which is signif- 
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FIG. 5: Revisited analysis of the historical s(Q,t) data on 
CuMn 5% (5). Below T 9 (~27.5 K) the NSE spectra were 
combined with values calculated from the dynamic suscepti- 
bility (10 4 > t > 10 7 ns). The data are plotted in a log-log 
scale and the continuous lines correspond to a simple power 
law decay below T 9 and to the Ogielski function above T s . 



icantly lower than % 2 = 0.55, the value obtained for the 
stretched exponential and \ 2 = 1-67 for the simple expo- 
nential respectively. The power law part of the Ogielski 
function, which holds at short times, describes the main 
part of the relaxation above T g and the spectra of Figure 
1 lead to an accurate determination of the effective power 
law exponent x as function of T . On the other hand, the 
stretched exponential only influences the tail of the relax- 
ation so the parameters t(T) and (3 are obtained here to 
low accuracy: r(T)« 0.3 ns and /3~1 at 55.8 and 50.8 K 
whereas t(T)r^22 ns and /3«0.66 at 45.7 K. In similar sys- 
tems, muon spin depolarisation measurements which are 
sensitive to the time range from 10 ns up to about 50 /xs 
(i.e. in the range where the relaxation above T g is mainly 
described by the stretched exponential) showed that in 
fact (3 approaches about 1/3 at T ^ 24 i 2 ^ as expected by 
the simulations. Susceptibility measurements very close 
to T g in an Ising spin glass show a similar limiting value of 
Furthermore, recent magnetic-field dependent muon 
spin depolarisation measurements on several AgMn spin 
glasses, analysed assuming an Ogielski-like decay of the 
correlations, lead to values of x tending to ~0.15 at T g 
in agreement with our results^ 

The large dynamic range covered by our data and 
their accuracy allow us to distinguish between a simple 
stretched exponential decay, the Ogielski function with 
(3 = 1 (power law times simple exponential), and the full 
Ogielski function; the full function is needed to give an 
acceptable fit of the data over the whole time and tem- 
perature range. 

In spin glasses the normalized s(Q,t) does not vary 
with Q, in dramatic contrast to the strong Q-dependence 
of the dynamics in ferromagnets. This remarkable Q- 
independence of the relaxation in spin glasses^ implies 



that the NSE s(Q,t) can be identified with q(t) and anal- 
ysed in the frame of dynamic scaling. The NSE spec- 
tra of Aun.86Feo.i4 collected on IN15 for Q — 0.4 and 
0.8 nm" 1 and on SPAN for 0.6 nm" 1 <Q<2.6 nm^ 1 con- 
firmed this behaviour. Fig. 3 compares, in a lin-log plot, 
spectra collected at 45.6 and 40.6 K for Q= 0.04 A^ 1 and 
0.08 A -1 respectively. The data almost overlap although 
the Q-values differ by a factor of two and the magnetic 
intensity decreases by almost a factor of 3, which ex- 
plains the larger error bars of the data set at 0.08 A -1 . 
The Q-independence of s(Q,t) was an important check of 
the high quality of our samples in a concentration range 
close to ferromagnetism. It also implies that s(Q,t) can 
be directly related to the macroscopic ac susceptibility 
X(Q=0, ui) = S(Q=0)[l-s(Q=0,t)]/kT,±iThe meaning of 
this equation is quite simple: s(Q, t) is the fraction of the 
total magnetic response S(Q) which does not relax be- 
fore the time t, in other words s(Q, t) is the part of S(Q) 
which cannot respond to a driving field of frequency 1/t. 
Figure 4 shows the NSE spectra for T < T g combined 
with s(Q,t) values deduced from macroscopic dynamic 
(a.c.) susceptibility measurements on the same sample. 
The data follow the power law decay over an impressively 
large range of at least 9 orders of magnitude in time. 

The impressive similarity between the experimental 
s(Q,t) and the decay of q(t) found in large scale simula- 
tions and more particularly the simple power law decay 
found at and below T g incited us to revisit the historical 
CuMn 5% data. Below T g , the spectra of CuMn 5% were 
combined with macroscopic dynamic (ac) susceptibility 
measurements and thus covered more than 9 orders of 
magnitude in time down to the microscopic time scaled 
These data, plotted now on a log-log scale, also reveal 
a power law decay below T g over an impressively large 
range of at least 9 orders of magnitude in time as shown 
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FIG. 6: Temperature dependence of the exponent x obtained 
by fitting the NSE spectra to a simple power low decay below 
T 9 and to the Ogielski function above T g . 
predicted by dynamic scaling. 
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in Fig. 5. Above T g , which is at 27.5 K, the power law 
holds only at short times and the decay can be described 
by the Ogielski function with /3«1 just as in AuFe 14%. 

The values of the exponent x, plotted versus the re- 
duced temperature T/T g in Fig. 6, are similar for both 
metallic systems. This confirms the similarity of the 
dynamic behaviour of AuFe, AgMn and CuMn systems 
around T s , pointed out by Uemura after comparing the 
NSE spectra of CuMn 5% with ^SR spectra of AgMn 
and AuFei^i As already mentioned, dynamic scaling re- 
lates the value of x at T g to the static and dynamic crit- 
ical exponents, which can be determined completely in- 
dependently from non-linear and ac macroscopic suscep- 
tibility measurements. Well established values of these 
exponents are available for a AgMn spin glass 6 and for 
CuMn and AgMn spin glasses doped with Auj£ in agree- 
ment with more recent data on AuFe 14%. The static 
exponents lead to rj— 2-j/is = 0.23±0.3 and the dy- 
namic exponent z=5.3±0.8 was determined by low fre- 
quency susceptibility measurements. From these val- 
ues we calculate x=0.116±0.026 in excellent agreement 
with x=0.116±0.007, the value we have observed on 
Auo.86Feo.14 at T g . It is important to note that the 
value of x in these basically Heisenberg spin glasses 
is considerably larger than that seen numerically^ or 
experimentally 2 ! in Ising spin glasses where x«0.07. 

Our results represent the first verification of scaling in 
spin glasses relating quantities of very different nature 
and measured by different methods on very different 
time scales, namely the microscopic time dependence 
of the autocorrelation function at times between the 
microscopic time and 10 -7 seconds on the one hand, 
and the macroscopic a.c. susceptibility at time scales 
greater than 10~ 3 seconds together with the static 
non-linear susceptibility on the other. This agreement 
constitutes a most compelling evidence for a true phase 
transition in spin glasses at T g . In addition, this demon- 
stration of critical relaxation characterisic of a phase 
transition also implies that the entire non-exponential 
temporal relaxation in these spin glasses is an intrinsic, 
homogeneous feature, as already evidenced by Uemura 
after comparing /xsR, NSE and dynamic susceptibility 
datai2& This conclusion is in obvious contrast to the 
clear evidence for the heterogeneous origin of the 
non-exponential relaxation in some structural glasses, 
where also no clear signature of a phase transition could 
be founds The clear evidence for a phase transition 
with a non-conventional order parameter and intrinsic, 
homogeneous non-exponential relaxation in spin glasses, 
a member of the large family of glassy systems, which 
also include structural glasses, other non-crystalline 
systems and living matter, is of particular importance 
in view of the unsolved central question of the nature of 
the glass transition in general. 
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